Tamm surface resonances in very low energy electron scattering from clean metal 

surfaces 
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Very low-energy features which occur in electron reflectivities from clean fee metal (111) surfaces 
have been subject to a number of interpretations. Here we analyse the feature near 19.8 eV on 
Cu(lll) at normal incidence and find that it is due to resonant scattering at the rise of the muffin- 
tin average interstitial potential between atomic layers on approach to the surface from the bulk. 
This new mechanism corresponds to a Tamm-type surface resonance which is very different in 
formation to the usual Shockley and Rydberg resonances and explains all features in a systematic 
way. 

PACS numbers: 61.14.Hg, 68.37.Nq, 73.20.At 



The copper (111) surface is of special interest because 
it has a band gap between the Fermi level and vacuum 
level where a Shockley surface state exists. Adsorbed al- 
kali metals on this surface produce interface states which 
also lie in this region. These states form a quasi 2D 
free-electron gas within a quantum well and it has been 
suggested that such systems can be used for device ap- 
plication at room temperature 0. It is therefore impor- 
tant that properties of the clean metal surface are com- 
pletely determined. These can be studied by photoemis- 
sion (PE), inverse photoemission (IPE) spectroscopies 
and other techniques which give (mostly) below vacuum 
energy level data. Above-vacuum energy-level data can 
be found from very low energy electron scattering spec- 
troscopies such as very low energy electron diffraction 
(VLEED), low energy electron microscopy (LEEM) and 
target current spectroscopy (TCS). Both regions can be 
analyzed theoretically by a layer-by-layer KKR scatter- 
ing method Q and properties determined from both en- 
ergy regions must exhibit a smooth continuation over the 
whole energy range. This approach effectively provides 
an expanded data base from which to unravel the compli- 
cated surface properties, many of which are expected to 
have significant energy and momentum variation. Infor- 
mation obtained from the higher energy region also has 
relevance for analyzing the below-vacuum region. 

Here we analyze VLEED and LEEM data on Cu(lll). 
The experimental data at normal incidence has a peak 
at ~ 20 eV which has not been easy to explain. It has 
variously been attributed to a band structure effect Q , a 
sub-threshold effect |4| and anisotropic inelastic scatter- 
ing 0. Similar features are also present for other (111) 
metal surfaces VLEED and LEEM data may have 
been under-utilized because of the difficulty in systemat- 
ically accounting for all the features. This is unfortunate 
because it is very sensitive to surface scattering poten- 
tials and electron self-energy including their energy and 
momentum variation. It is also more sensitive to the 
vertical position of surface atoms than LEED which is 



performed at energies higher than ~ 40 eV. 

Shockley and Rydberg states and resonances have been 
detected in below-vacuum PE and IPE experimental data 
and resonances in above- vacuum VLEED, LEEM and 
TCS data. Historically surface states (and resonances) 
have been labeled as Tamm-type and Shockley-type al- 
though in some contexts there is no sharp distinction. In 
a scattering approach, the surface state can be labeled a 
Shockley type if an essential element in its formation is 
the shape of the surface barrier in the vacuum beyond the 
top row of atoms. (They have also been called crystal- 
induced surface barrier states in this context). Rydberg 
states arise from the image potential tail of the surface 
barrier beyond the top row of atoms and are also called 
image surface barrier states. In this context a Tamm 
surface state/resonance could be defined as one which 
arises from change in the 3D periodicity of the crys- 
tal near the surface but is not of the Shockley type as 
defined above and not due to adsorbed foreign surface 
atoms. In a scattering picture, the Shockley and Ryd- 
berg states/resonances arise because of standing waves 
forming at some energies in the potential well consist- 
ing of the crystal on one side and the rise of the surface 
barrier to the vacuum on the other. 

With the vacuum emergence of the {01}, {10} beams 
at 30.7 eV for Cu (111) at normal incidence and an ex- 
pected crystal inner potential of ~ 13 eV .4], the 20 eV 
peak lies in the energy range between the crystal and vac- 
uum emergence of the above non-specular beams. In this 
case scattering of these beams at the surface barrier can 
produce Shockley resonances. If this situation were the 
case, the 20 eV peak would give important information 
about the form of the surface barrier near where it joins 
to the crystal which has consequences for the determina- 
tion of other barrier properties. Such a low-energy lying 
Shockley-type resonance has not been identified to date. 

This possibility has been examined but the 20 eV peak 
width, height and Bragg-peak separation could not be 
reproduced by this mechanism for realistic values of elec- 
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FIG. 1: Calculated reflectivity of 00 beam on Cu(lll) at nor- 
mal incidence with inner potential U Q , uniform inelastic ab- 
sorption potential U ln , crystal temperature Tand Debye tem- 
perature 315 K. Frames (a), (b), (c) use a non-reflecting sur- 
face barrier. Frame (d) has an image barrier at z = —5.5 a.u. 
with saturation to 0.4 eV at —4.9 a.u. and barrier inelastic 
absorption described in Ref. 8 with a = 1.33 a.u. The down- 
ward/upward arrows indicate vacuum/crystal emergence of 
the {01}, {10} beams respectively. 



tron inelastic scattering of ~ 2.5 eV |4],|6j. One model 
came close with the resonance forming a dip in the Bragg 
peak splitting it into two peaks. This is shown in Fig. 1. 
However the position of the image origin had to be placed 
at distances greater than that which is expected theoret- 
ically. At realistic positions of the image origin, Zq> (< 4 
a.u. from the centre of the top row of atoms), the well was 
not wide enough to support a strong resonance only a few 
eV above its base for electron absorption of U ln — ~ 2.5 
eV. The bulk elastic scattering potential used in all the 
calculations here is due to Moruzzi et al 0, and it was 
found to give excellent agreement with band gaps below 
the vacuum level. 

Another means by which a resonance could occur a few 
eV above the crystal-emergence energy of the {01}, {10} 
beams is because of a rise in the muffin-tin average in- 
terstitial potential between atomic layers from the bulk 
value on approach to the surface. This is illustrated in 
Fig. 2. Such changes must occur because of the change 
in environment of surface atoms. No features due to this 
change in the 3D periodicity of the crystal appear to have 
been identified in low energy electron spectroscopies to 
date. From a scattering point of view, a potential well 
is set up where standing waves can be formed between 
the crystal Bragg pseudo-gap (strong bulk reflection) and 




Distance from centre of top row of atoms z (a.u.) 

FIG. 2: Schematic potential energy at the Cu(lll) surface 
in a direction perpendicular to the surface. Symbols are ex- 
plained in the text except for Zj, the jellium discontinuity. 
Dashed horizontal lines with double-ended arrows near Z\ in- 
dicate electron scattering which gives rise to standing waves 
and the feature near 20 eV in experimental data. Unshaded 
regions in the substrate represent surface-projected bulk-band 
gaps. 



the potential rise. In a scattering picture, this is how a 
Tamm-type surface resonance is formed. Such a rise in 
interstitial potential could also produce Tamm surface 
states below the vacuum level in surface-projected bulk 
band gaps and resonances outside these gaps. The point 
is that potential variations detected in above-vacuum 
spectroscopies (like VLEED, LEEM) also have analogies 
for below- vacuum spectroscopies (PE, IPE etc). 

To examine this case, calculations were performed with 
an interlayer potential rise of height U\ — 0.1029 Ry = 
1.4 eV with respect to the bulk muffin-tin zero (crystal 
inner potential) placed between the top layer of atoms 
and the atomic layer below. This is illustrated schemat- 
ically in Fig. 2. The top layer was followed by a surface 
barrier of height U-x = 9.75 eV with respect to U\. It 
was found that the shape of the surface barrier was not 
important here since it does not give rise to any barrier 
features of the Shockley or Rydberg type when the image 
origin is placed at any realistic position from the centre 
of the top layer of atoms and realistic crystal electron 
absorption of > 2.3 eV is used. The layer KKR method 
for calculating LEED reflectivities of McRae @ was used 
where the crystal substrate is terminated at z = and z 
is the distance perpendicular to the surface and positive 
directed into the crystal. Cu(lll) has a lattice constant 
of 6.8309 a.u. with 3 layers in the repeating unit. The 
top layer of atoms is taken as part of the selvedge and 
for Cu(lll) is placed at z = —3.9438 a.u. and trans- 
lated with respect to the first substrate layer which has 
an atom at the origin by x = 1.6102 and y = 3.2201 
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FIG. 3: Reflectivity of 00 beam on Cu(lll) at normal incidence. Frame (a) is the experimental result of Ref. 5 with energy 
step 0.5 eV. Frames (b) — > (h) are calculated profiles with U = 11.15 eV and Ui n = 2.5 eV in all layers except frame (f) 
which has Uin = 0.5 eV in the top atomic layer. Symbols are described in Fig. 1. In all cases an interlayer step of height 
U\ — 1.4 eV is placed between the top layer and the next layer at positions z\ from the centre of the second row of atoms at 
the surface, as indicated on the diagrams. The surface barrier height is 9.75 eV, the crystal temperature is 300 K and Debye 
temperature 315 K. Downward/upward arrows are the same as in Fig. 1. Dashed upward arrows indicate crystal emergence 
energy of {01}, {10} beams in the top atomic layer. Frames (g) and (h) have inelastic absorption at the interlayer step as 
described in the text. 



a.u. The interlayer rise in potential is modeled in these 
calculations by a simple step potential placed between 
z = and -3.9438 a.u. From higher energy LEED anal- 
yses on Cu(lll) the surface atomic layers have not been 
found to deviated significantly from their bulk positions 
4J. From the position z — —3.9438 a.u. there is a further 
rise in potential representing the surface barrier of height 
U2- The total rise in potential from the bulk muffin-tin 
zero to vacuum level is U\ + U2 = Uq = 11.15 eV in the 
calculations shown here. 

Wave vectors of the wave functions corresponding to 
beam v in the three regions bulk, selvedge and vacuum 
are represented by fe*, kf and respectively after 
the notation of McRae Jj|. For the interlayer step, the 
reflection and transmission coefficients are 
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corresponding wave vectors and z\ is the position (with 
respect to the origin) of the potential step. This sharp 
step overestimates above-step reflections but the aim is 
to keep the model as simple as possible at this stage and 
easily reproducible by others. The ps and t$ were used 
to calculate a transfer matrix Q2 for this step scatter- 
ing layer according to Eq. (14) of McRae [9j. Similarly 
transfer matrices for the top atomic selvedge layer Qi 
and for the surface barrier layer Qo were calculated, giv- 
ing a selvedge scattering matrix X where 



X — Q2 ■ Qi ■ Qq. 



(2) 



where k^~, 



are the perpendicular components of the 



The selvedge scattering matrix X was combined with 
the bulk crystal scattering matrix M using Eq. (35) of 
McRae . In order to keep the model as uncomplicated 
as possible, the phase shifts for the top atomic selvedge 
layer were kept the same as the bulk. The bulk isotropic 
inelastic scattering potential [/i n (buTk) was 2.5 eV and 
for the selvedge atom layer /7i n (selv) was chosen from 0.5 
to 2.5 eV. Different heights U\ of the interlayer step (and 
barrier height U2) were tried. 

Fig. 3 shows the results of our calculation for inter- 
layer step positions from Z\ = -0.7 to -3.7 a.u. with U\ — 
1.4 eV. In frames (b) to (d), no inelastic scattering was 
included for the step layer, i.e. the normal components 
of the wave vectors in Eq. (1) had no imaginary compo- 
nents. This of course overestimates the effects but allows 
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one to see the possibilities. Firstly, the value of inelastic 
scattering in the selvedge top layer, {/^(selv), has some 
effect on the fine structure features but for this simple 
model, optimization of such parameters was not consid- 
ered appropriate at this stage. In comparison with the 
reflectivity without the interlayer step as seen in frame 
(a) of Fig. 1, we see that a peak appears near 20 eV 
and moves down in energy from 20.3 to 20 eV as \zi\ 
increases. The six {01}, {10} beams emerge in the bulk 
crystal layers at 19.5 eV and in the selvedge top atomic 
layer at 20.9 eV (with respect to the vacuum level). The 
{01}, {10} beams from the bulk are propagating and in- 
cident on the step in a direction towards the surface in 
the energy range above 19.5 eV. The {01}, {10} beams 
in the top atomic layer are incident on the step in a di- 
rection away from the surface and are evanescent until 
20.9 eV when they become propagating. Multiple reflec- 
tions and transmissions in both directions give rise to the 
{01}, {10} beams becoming incident on the bulk atomic 
layers and selvedge atomic layers with phases and am- 
plitudes different from what occurs without the step. At 
certain energies, standing waves may form which corre- 
spond, in a scattering picture, to what historically has 
been called Tamm surface resonances. The non-specular 
scattering at atomic layers leaks some flux into the 00 
beam. Above step-height reflections may be exaggerated 
here but some will also occur for a smoother interlayer 
potential rise. 

In this model, inelastic scattering in the whole volume 
of the selvedge and substrate regions already has been 
taken into account in the calculation of the bulk and 
atomic layer scattering matrices. Also inelastic scatter- 
ing in the top bulk layer is likely to be different from 
that for lower positioned layers. It is difficult for a model 
based on the analytical reflections coefficients in Eq. (1) 
to ascribe appropriate inelastic scattering from events as- 
sociated with the step only. Changes in the isotropic- 
layer inelastic scattering potentials on approach to the 
surface are unknown quantities at this stage. Hence for 
this model we have reduced the intensity of the interlayer 
step features somewhat arbitrarily through the wave vec- 
tors in Eq. (1) by adding an inelastic scattering potential, 
C/j n (step), of 0.3 eV and 0.15 eV. This simulates how the 
features change with inelastic scattering. We see in frame 
(e) of Fig. 3 that the structure appears to be of the same 
form as seen in the experimental data including its energy 
width and separation from the Bragg peak. Resolution 
of any fine structure features would also be lost in most 
experimental set-ups because of beam width and energy 
broadening. With exaggerated above-step features the 
model used here would not be expected to exactly repro- 
duce the experimental result, but the origin of the feature 
can be determined. 



In conclusion we have found the following: (a) For nor- 
mal incidence electron reflectivity from Cu(lll), a peak 
near 20 eV with the given experimental width, height 
and Bragg-peak separation could not be produced from 
a surface barrier with image origin, zq, located in a real- 
istic position and for inelastic scattering potential, Ui n , 
of realistic value; the peak is not due to a Shockley-type 
surface barrier resonance, (b) The 20 eV peak is due to 
resonant scattering at the rise in the average interstitial 
potential between atomic layers on approach to the sur- 
face. It is therefore a Tamm-type surface resonance. This 
type of surface resonance has not been identified before 
in VLEED, LEEM etc. It is also likely to account for 
similar features in other clean metal (111) surfaces, (c) 
The interlayer rise in potential may only produce surface 
resonances when they fall in the energy range of a surface- 
projected bulk band pseudo-gap (Bragg peak) in VLEED 
and related spectroscopies, (d) Even when the condition 
in (c) above is not fulfilled, the rise in potential before 
the top atomic layer will change the value of the height of 
the surface barrier. This value affects the determination 
of barrier features such as image origin position, satu- 
ration etc, for both below vacuum level spectroscopies 
(PE, IPE) and those above (VLEED, LEEM, TCS). (e) 
All of the features in VLEED (0 - 40 eV) may now be 
accounted for in a systematic way if, (i) scattering at a 
realistic image surface barrier, and (ii) scattering at the 
interlayer rise in interstitial muffin-tin potential between 
surface layers is included in the theoretical calculation of 
reflectivities. This would make VLEED and LEEM more 
widely used tools for unravelling surface properties in the 
future. 
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